Unified Description of Nuclear Stopping in Central Heavy-ion Collision from 10A 

MeV to 1.2 A GeV 
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The detailed analysis of wide excitation function of nuclear stopping has been studied within a 
transport model, Isospin-dependent Quantum Molecular Dynamics model (IQMD) and an overall 
good agreement with the INDRA and FOPI experimental data has been achieved. It is found that 
mean value of isotropy in central Heavy-Ion Collision (HIC) reaches a minimum near Fermi energy 
and approaches a maximum around 300 - 400A MeV. This suggests that, in statistical average, the 
equilibration is far from being reached even in central HIC especially near Fermi energy. A hierarchy 
in stopping of fragments, which favors heavy fragments to penetrate, provides a robust restriction 
on the global trend of nuclear stopping and could serve as a probe for nuclear equation of state. 

PACS numbers: 25.70.-z, 21.65.Mn 
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I. INRODUCTION 



The equation of state (EOS) of nuclear matter and its 
transport mechanisms are the focus of Heavy-ion Col- 
lisions (HIC) during the past three decades [lrla]. Be- 
fore the early 1980s, hydrodynamics approaches based on 
the local equilibration postulate were provided to extract 
the EOS information from experimental observables Q. 
The discovery of collective flow seems to confirm these 
approaches and their corresponding equilibration condi- 
tions. With the development of microscopic transport 
theories, it is found that the local equilibration postulate 
is not indispensable. Actually, near Fermi energy, both 
statistical and dynamical models which are based on very 
different, even contradictory mechanisms can predict the 
same multifragmentation phenomenon (lfjj . 

During central HIC process, the nuclear stopping gov- 
erns most of the dissipated energy and constrains the 
different reaction mechanisms at different incident ener- 
gies. It can provide the information on the EOS, nucleon- 
nucleon (N-N) cross section and the degree of the equi- 
librium reached in HIC {[MM- Recently, INDRA and 
ALADIN Collaborations investigated event- by- event nu- 
clear stopping in central HIC at intermediate energies by 
4-7r multidetector [12]. In the experiment [r|, a striking 
minimal stopping at AOA MeV (Fermi energy) for Xe+Sn 
system was observed. These experimental results are po- 
tential to provide important information about the prop- 
erties of the nuclear matter and help to find out critical 
clues for the whole dynamical process during HIC. In pre- 
vious comparison between experimental results and one 
isospin-dependent Quantum Molecular Dynamics model 
(IQMD) simulations [l3| , a significant difference between 
them is observed. In the experiment [HJ, below 70 A 
MeV, they show the experimental data far below the 



IQMD simulation results and argued that IQMD model 
overpredicts nuclear stopping power at low energy. How- 
ever, before comparing the experimental results quan- 
titatively with that of transport model, two important 
ingredients should be carefully considered. One is the 
difference between nucleon phase space (the momentum 
and positions ensemble of all nucleons) and fragment 
phase space (the momentum and positions ensemble of 
all fragments) in IQMD, and the other is the criterion of 
the centrality in the experiment. After considering these 
two gredients, we find that the nuclear stopping could be 
served as a potential probe for determining the nuclear 
EOS. 

This article is organized as following: After introduc- 
ing some details on IQMD model in Sec. II, two gredients 
are reconsidered in Sec. Ill: one is the difference between 
nucleon phase space and fragment phase space and the 
other is the impact parameter mixing. Then we com- 
pare the simulation results of IQMD with INDRA's ex- 
periments near Fermi energy as well as FOPI's at higher 
energy in Sec. IV. Summary and conclusions are given in 
Sec.V. 
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II. ISOSPIN DEPENDENT QUANTUM 
MOLECULAR DYNAMICS (IQMD) MODEL 

The Quantum Molecular Dynamics (QMD) [2l[ model, 
is a n-body theory, to predict the behavior of HIC at 
intermediate energies on an event by event basis. The 
IQMD [22|,[23[ model, as inheriting from the QMD model, 
considers the isospin effects at various aspects: different 
density distribution for neutron and proton, the asym- 
metry potential term in mean field, experimental cross- 
section for necleon-necleon (a np « 3<7 pp ,(7pp = a nn ) and 
Pauli-blocking for neutron and proton separately. Just 
like in QMD, each nucleon wave function is represented 
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as a Gaussian form in IQMD, which defined as: 
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Here L is the width parameter for the Gaussian wave- 
packet, which is found to be dependent on the size of the 
reacting system, to constrain the stability of system. Wc 
use L = 2.16 fm 2 for Xe+Sn system. The fl(t) and Pi(t) 
are the position and momentum coordinates of nucleon. 
Performing variation method, equations of the time evo- 
lution of the mean position rl{i) and momentum pi (t) are 
found to be well known Hamilton equations of motion: 



Pi 



d(H) 



and 



d(H) 



(2) 



where (H) is the Hamiltonian of the system. So now the 
problem is to calculate the Hamiltonian of the system. 

After applying the Wigner transformation to the single 
nucleon wave function, one can get Wigner distribution 
fuction to describe a single nucleon density in phase space 
defined as: 



fi(r,p,t) 
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and the total density of the system is the sum over all 
the nucleons. 

Then we can get total Hamiltonian in IQMD model: 



(H) = (T) + (V), 



(4) 



where the meanfield part is 



W = \ E E / M p >& f ) V ij fj(r',P', t) drdr'dpdp 1 . 

(5) 

In QMD model, the two-body potential interaction con- 
tains the Coulomb interaction and the real part of the 
G-Matrix. The later one can be divided into three 
parts: the contact Skyrmc-type interaction, a finite range 
Yukawa-potential, and a momentum dependent part. 
Mathematecally, the two-body potentail interaction can 
be written as: 
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where Z is the charge of the baryon, t\...ts and p are 
the parameters to fit the real part of G-matrix and the 
properties of nuclei. For IQMD model, one additional 
part, the symmetry potential V^ m , is added to take into 



account the isospin effects. The total potential in IQMD 
is then written as: 

V — ^Skyrme ^ ^Yuk + " mdi ^ V Cou\ ^ v symi \ ' ) 

where the symmetry potential 



Kym = t(i-^T 3i T 3j S(r t - rj) 



(8) 



with = 100 McV for fitting the Bethc-Wcizsacker mass 
formula, T3 is the isospin third projection and po is the 
nuclear saturation density (0.16/m~ 3 ). 

We focus on the Skyrme potential, the momemtum 
dependent potential, due to their special connection to 
nuclear EOS. 

One can fulfill the calculation for Skyrme potential and 
momemtum dependent in Eq. [5] by introducing the in- 
teraction density, 



Pn 



The momentum dependent potential, which may be 
optional in QMD and IQMD, is parameterized with the 
experimental data, can be written as 
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where 5 = 1.57 MeV and e = 500 (GeV/c) 2 are taken 
from the measured ener gy dependence of the proton- 
nucleus optical potential |2J] and pint — X)/°int( 7 ^)- 
The Skyrme potential is 



{/Sky = a( Wntj + p^Pint ^ 

Pa Po 



(11) 



where a, (3 and 7 are the Skyrme parameters, which con- 
nect tightly with the EOS of bulk nuclear matter. After 
fitting the minimum binding energy and the compressibil- 
ity at the saturation density, one can get the parameters. 
The nuclear compressibility, which is the second deriva- 
tive of the energy at the minimum with respect to the 
density is expressed as: 



i:) 2 



E 



(12) 



K = 200 McV means soft EOS, while k = 380 McV is for 
hard EOS. Please see Table U the Skymre sets of param- 
eters for different EOS. 

The other terms can be also calculated by applying 
the convolution of the interaction with the Wigner dis- 
tribution function. After all the potential terms are de- 
termined, one can solve the equations of motion of the 
baryons. For the mesons, only Coulomb force is consid- 
ered. 

For the collisions, IQMD uses the experimental cross- 
section which tell the isospin effects and nuclear medium 
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a (MeV) (MeV) 


7 


8 (MeV) e 


(cevO 


s 


-356 303 


1.17 






SM 


-390 320 


1.14 


1.57 


500 


H 


-124 71 


2.00 






HM 


-130 59 


2.09 


1.57 


500 



TABLE I: Parameter sets for the nuclear equation of state 
used in the QMD model. S and H refer to the soft (compress- 
ibility k —200 MeV) and hard equations of state (compress- 
ibility k =380 MeV), M refers to the inclusion of momentum 
dependent interaction. The results are taken from [S^] 



effect [22|. To consider the fermion property, the Pauli- 
blocking is also adopted after the collisions. 

IQMD treats the many body state explicitly, contains 
correlation effects to all orders and deals with fragmen- 
tation and fluctuation of HICs natually. To recognize 
fragments produced in HICs, a simple coalescence rule 
is used with the criteria Ar = 3.5 fin and Ap = 300 
MeV/c between two considered nucleons. Thus, nucle- 
ons dominated in Fermi motion will be constrained in 
the fragments. 
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FIG. 1: (Color online) Comparison between experimental stopping 
values [12| | (black points) and predictions of IQMD calculation in 
nucleon phase space (red open circle) and fragment phase space 
(blue up triangle) for central 129 Xe + 120 Sn collsion. Two lines are 
for guiding the eyes. IQMD adopts the hard EOS and fe ret ; < 0.1. 



III. DIFFERENCE BETWEEN NUCLEON 
PHASE SPACE AND FRAGMENT PHASE 
SPACE AND IMPACT PARAMETER MIXING 

To describe nuclear stopping power, the ratio of trans- 
verse to parallel quantities is used in the experiment fi~2| . 
One definition is the energy-based isotropy ratio Re, 



Re = 



(13) 



where Eti(Eu) is the cm. transverse (parallel) kinetic 
energy and the sum runs over all products event by event. 
One can expect Re = 1 for isotropy or full stopping, 
Re < 1 for partial transparency, and Re > 1 for super 
stopping. Another definition is R p , 
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where momenta are used instead of energies. However, 
they are actually the different forms of physics realization 
of classical Maxwell distribution assumption. 

At low energy, two-body collision between nucleons is 
suppressed by Pauli-blocking greatly and the mean field 
governs the HIC process. At freeze-out stage, the motion 
of a large portion of nucleons is limited in certain frag- 
ments, in which Fermi motion dominates over the others. 
While, at higher energy, two-body collision comes to play 
the dominant role and most of nucleons suffer violent col- 
lisions and get excited to become free. In this case, the 
fragment phase space tends to get closer to the nucleon 
phase space. While, in previous IQMD simulation by Liu 
et al. [l3j, the nuclear stopping was investigated from 
nucleon phase space at intermediate energy, where the 



isotropic Fermi motion still plays a great role. In com- 
parison with the experimental data of the stopping, the 
fragment phase space should be applied at intermediate 
energy within IQMD model rather than nucleon phase 
space. 

In Fig. [TJ we show the results of stopping defined as the 
ratio of transverse momenta to parallel momenta R p , in 
nucleon phase space and fragment phase space, respec- 
tively. The scaled impact parameter, b rec i < 0.1 (b rec i = 



is adopted, with b n 



1.12(Ap /3 



A]{ 3 ) fin, 



where Ap and At arc the mass of projectile and target, 
respectively. Apparently, stopping calculated in nucleon 
phase space reaches almost the saturated value, just like 
Liu's situation [l3j . While, in the case of fragment phase 
space, stopping become far below the saturated value and 
get closer to the experimental value. It underlies that, 
at low incident energy, nuclear stopping obttained from 
fragment phase space is more sensitive as compared to 
the nucleon phase space and tends to fill the gap between 
theory and expeiments. At higher incident energy (above 
100 A MeV), the difference between stopping presented in 
nucleon phase space and fragment space in IQMD, tends 
to get smaller. However, a minimum stopping with value 
0.70 at about 55 A MeV is shown in IQMD, compared 
that with value 0.68 at A0A MeV in experiment. This 
difference may come from the different criteria of cen- 
trality between IQMD and experiment. 

To avoid autocorrelations, total charged particle mul- 
tiplicity (N c h) in the experiment [l2[ is chosen as the 
criterion of the centrality. However, the largest N c h does 
not mean the highest centrality, due to event- by- event 
fluctuation. On one hand, N c h gets saturated at certain 
centrality and then central and near central HIC provide 
similar N c h- On the other hand, due to the dynamical 
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FIG. 2: (Color online) IQMD simulation with the hard EOS for 
129 Xe + 119 Sn collisions at 50A MeV. (a) Bi-dimensional cor- 
relation plot between the mean value of stopping and the total 
charged particle multiplicity, same as Fig. 1(b) in Ref. Q. (b) Bi- 
dimensional correlation plot between the mean value of stopping 
and impact parameter, (c) Distribution of Re for central colli- 
sions, (d) same as (b), but with central collision cut. 
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FIG. 3: (Color online) Energy-based stopping values evolution 
with time for 129 Xe + 120 Sn at 50A MeV and b red < 0.15 in IQMD 
(open symbols) and experimental data (solid line) [T3 |. Circles, 
squares, up-trianglcs and down-triangles represent our calculations 
with the hard (H), the soft (S), the hard with MDI (HM), and 
the soft with MDI (SM) for EOS, respectively. The horizontal line 
represents the measured result by the INDRA collaboration (not 
for the time evolution). 



IV. RESTULS AND DISCUSSIONS 



fluctuation, N c h at certain centrality presents a broad 
distribution. Thus, the HIC at different ccntralities may 
share the same N c h with a great chance. Hence, if N c h is 
adopted for the criterion of the centrality, it is expected 
that the impact parameters will suffer a great mixing 
and cover a wide range, which will decrease the stopping 
power and increase the fluctuation of the HIC system 
greatly. 

We now use IQMD model to elucidate the impact pa- 
rameter mixing situation. In the experiment [12j, a typ- 
ical cross-section of 50 mb for the most total charged 
particle multiplicity is taken for all HIC systems. Within 
the IQMD model, we calculate the mean N c h and its 
width at head-on collision (b re d = 0). Then, low limit is 
set as the mean N c h with its half width off. In Fig. [2£a), 
IQMD model reproduces well the experimental results of 
INDRA and ALADIN Collaborations Q3], except giving 
a little higher N c h (note that the experimental filter is 
not used in the present calculation). The red dash line 
at N c h = 44 is the low limit for central collision in IQMD 
which is comparable with N c h = 36 in experiment. In 
Fig. (lib), a mean value of stopping (Re) around 0.6-0.7 
is observed in the exact central collision, which has a 
decreasing trend with the increase of impact parameter. 
When the cut N c h — 44 is applied, there is a broad dis- 
tribution of impact parameter in Fig. [2Id), and a lower 
mean value of stopping about 0.56 with a width about 
0.42 in Fig. [5Jc), which is comparable with that a mean 
value 0.56 with a width 0.47 in experiment (T^]. In this 
context, better criterion of centrality should be provided 
to reduce the probability of impact parameter mixing at 
low energy. 



Taking the above two key ingredients into account, 
we can compare the experimental results with those of 
IQMD quantitatively. To have a broad view of the nu- 
clear stopping, a systematic simulation of 129 Xe + 120 Sn 
is made in very central HIC (b re d < 0.15) from 10^4 MeV 
to 1.2 A GeV to cover most of the intermediate energy 
range, where our results can be compared with the re- 
cent INDRA and FOPI data [HI, [T3]. 

It should be noted that related but different definitions 
of stopping are adopted by INDRA [U] and FOPI [H 
li~7l Il9j Collaborations, respectively. The event- by- event 
stopping (event level) defined as the ratio of transverse 
to parallel kinetic energies is adopted by INDRA and 
ALADIN Collaborations (i.e. Re)- While, the ratio of 
the variances of the transverse to that of longitudinal 
rapidity distribution of fragment (fragment level) is taken 
in the case of FOPI experiment 



vartl = 



(AyQ 2 
(Ay ; ) 2 ' 



(15) 



Despite this, they are the same in spirit to show the 
global stopping power of central HIC. When a full stop- 
ping stage (likely to be in equilibration) is reached, both 
Re and vartl are required with a value of unit. 

Now we compare our results with INDRA experiment. 
The central (b red < 0.15) collision of 129 Xc+ 120 Sn, at 
incident energy 50A MeV is taken as an example. The 
time evolution of the stopping (defined in Re) for vari- 
ous EOS are shown in Fig. [3l The stopping for different 
EOS presents the same trend. The stopping rises up at 
the early compression stage, drops down at the follow- 
ing expansion stage and becomes stable after the system 
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reaches the freeze-out. This is consistent with the results 
from works of [f| and 14 1, which are in nucleon phase 
space. However, the stopping power shows EOS depen- 
dence. The hard EOS shows higher stopping than the 
soft EOS, during the whole time evolution. This is due to 
the hard EOS nuclear matter is harder to be compressed 
than the soft one. When compressed from the longitudi- 
nal direction, the hard EOS nuclear matter will squeeze 
out more than the soft one, through the tranverse ex- 
pansion. The momentum dependent force decreases the 
stopping power of the system, this may result from its 
repusive nature and increasing in mean free path [2~H ]. 
At this incident energy, the INDRA experimental result 
seems to support the soft EOS. 

It may be pointed out that here the time adopted in 
IQMD for the freeze-out stage is 200fm/c. At low in- 
cident energy case, the heavy fragments may still stay 
at an exciting state, which would deexcite themselves 
by emissing free nucleons, LCP and gamma rays at the 
later stage. Basing on this consideration, some hybrid 
models, a dynamical one followed by a statistical one, 
come into the play, such as AMD + GEMINI [25ti27j and 
HIPSE + SIMON [Hi etc. They do fit the experimental 
fragment distribution very well. At the same time, the 
deexcitation procedure may also extend both transverse 
and longitudinal width of rapidity of the final fragments, 
which smears the original stopping at fragment level. Es- 
pecially, the big fragments at high excitation stage will 
also emit free nucleons or LCPs, which will increase the 
isotropy of the system. This might be very serious for 
lower energy (<40AMcV) case due to more heavy frag- 
ments are left. As a compensation, the freeze-out time 
is delayed from 200fm/c to 300fm/c (or even later), but 
no great changes happen, except for the IMF raise their 
stopping value a little higher. For higher energy (> 100/1 
MeV), since there exists very few heavy fragments, it is 
not necessary to evolute the system with such a long time 
or deexcite the system with a statistical process |29j . 

Fig. @] shows the excitation function of the mean value 
of Re and its width gr e for 129 Xc+ 120 Sn by using differ- 
ent EOS. In Fig. H[a), the soft EOS with the compress- 
ibility of k = 200 MeV is taken to show the distribution of 
Re with incident energies. In Fig.|4](b) and (c), the mean 
Re and its width ur b show the same trends, consistent 
with the experimental results [l2|. Scaled by 0.75, the 
width of stopping in the experiment of INDRA is com- 
pared with the simulation results. The broader widths in 
experiment could come from the impact parameter mix- 
ing, as we have discussed in Sec. II. For all the considered 
energy range, the hard EOS with k — 380 MeV shows 
stronger stopping than the soft EOS. The EOS with mo- 
mentum dependent interaction (MDI) [24[ tends to de- 
crease the stopping at energy below 500A MeV. While 
the weight of MDI tends to vanish above 500^4 MeV. 
Near the Fermi energy, we also observe a minimum of 
stopping with its minimal width for all EOS. The EOS 
with MDI favors more penetration and consumes less en- 
ergy to reach the minimal stopping. The present experi- 
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FIG. 4: (Color online) Stopping and its width as a function of 
incident energy with different EOS for central 129 Xe + 120 Sn colli- 
sion in IQMD (open symbols) and experimental data (solid circles) 
H2I . (a) Bi-dimensional plots show the distribution of Re- (b) 
The mean value of stopping, (c) The width of stopping (Note the 
widths of experiment have been scaled by a factor 0.75). Circles, 
squares, up-triangles and down-triangles represent our calculations 
with the hard (H), the soft (S), the hard with MDI (HM), and the 
soft with MDI (SM) for EOS, respectively. 



(b) 

s Ni+Ni b rpH <0.15 - 




NWS? 



V. A 



□ 



oH AHM 
h -«-lndra 

S, °S VSM 



(MeV/nucleon) 



10 



10 2 



10 J 



JMeV/nucleon) 



FIG. 5: (Color online) Same as Fig.g] but for 58 Ni + 5s Ni system. 
One more case (HM), is added in (c) to see the distribution of Kg. 



mental results seem to favor the soft EOS around Fermi 
energy. This also consists with recent conclusion on the 
soft hadronic matter basing on Kaon spectrum analysis 
at higher energy 0,0). At 300 - 400 A MeV, a maximum 
of stopping accompanied by its maximal width is seen for 
all EOS. This is comparable to the recent experimental 
results of FOPI collaboration [lj Q3, El]. 

The smaller size collision system 58 Ni + 58 Ni is also 
investigated to see the stopping behavior and similar be- 
haviors are found in comparison with the case of 129 Xe 
+ 120 Sn, except for 58 Ni + 58 Ni shows lower stopping 
power. It may be noticed again, the INDRA data seems 
to support the soft EOS case. (Fig. Ob)) 

Fig. [6]Ja) displays the atomic number (Z) hierarchy of 
the degree of stopping at fragment level at different inci- 
dent energies. From small fragment Z = 1 with stopping 
value 0.7 — 0.9 to intermediate one Z = 10 with 0.1 — 0.2, 
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FIG. 6: (Color online) Stopping hierarchy in IQMD with the 
soft EOS for central 129 Xe + 120 Sn collision. Four energy points, 
namely 40, 100, 300 and 1200,4 MeV are investigated, (a) and 
(c) shows Z-dependent stopping at fragment level and the yield of 
fragments (scaled by Z=l), respectively, (b) and (d) displays bi- 
dimensional distributions of transverse and longitudinal rapidity 
for different fragments at 40,4 MeV, respectively. 
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FIG. 7: (Color online) Atomic number dependence of the stopping 
observable vartl for fragments emitted in central 129 Xe + 120 Sn 
collisions. Open symbols represent our calculations with different 
EOS and solid symbols with error bar denote the available exper- 
imental data [Till , (a) 45,4 MeV (no experimental data available 
yet); (b) 150,4 MeV; (c) 250,4 MeV; (d) 400,4 MeV. 



indicates that the stopping power drops quickly as Z of 
the fragment increases. Similar trend has been experi- 
mentally observed Au + Au, Xe + Csl and Ni + Ni cen- 
tral collisions at 150 A MeVj_25(L4 MeV and 400 A MeV, 
by FOPI collaboration [O, [H . 

The hierarchies also show the incident energy depen- 
dence. On the one hand, at 300,4 MeV (open diamond 
in Fig. H][a)), fragments show stronger stopping than 
those at the others which is consistent with the maximum 
event-level stopping at the same energies in Fig. [4] This 
may suggest a state close to equilibration. On the other 
hand, at AO A MeV (open square in Fig. (B^a)) where the 
stopping at event level reaches a minimum (Fig. EJb)), 
the stopping at fragment level does show its correspond- 
ing smallest value. At this incident energy, most of nu- 
cleons are confined in certain fragments, which also re- 
duces the stopping value at event level. Near Fermi en- 
ergy, the relative yields of intermediate-mass fragments 
(scaled by Z = 1) arc higher than those at the others 
as shown in Fig. Htc), which will lead to the smallest ef- 
fective power-law exponents in charge distribution [3l| . 
For a further understanding the minimal stopping near 
Fermi energy, in Fig. EJb) and Fig. EJd), the holographic 
2D histogram of transverse and longitudinal rapidity for 
different charge of fragments are presented (rapidity has 
been scaled by the projectile rapidity in center of mass 
system). From Z = 1 to larger Z, the longitudinal rapid- 
ity distributions become relative broader than the trans- 
verse one. In addition, for those fragments with Z larger 
than about 4, two clear peaks in the longitudinal rapidity 
stand prominently, corresponding to projectile-like and 
target-like contribution. This means at Fermi energy, 
equilibration is far from being reached and entrance chan- 
nel effects remain strong even in central HIC. While, at 
high energy 1.2 A GeV, fragments are the most transpar- 
ent. 



It is highly expected that the hierarchy would help 
to determine nuclear EOS. Fig. [7] shows Z-dependent 
stopping for four different energies with different EOS 
in IQMD. The hard EOS case shows higher stopping at 
fragment level than the soft EOS one, and the MDI al- 
ways makes it easy for the fragments to penetrate. At 
lower energy (150 A MeV in FigJT^b)), the EOS without 
MDI agrees with FOPI data [llpmich better (Xe+Csl 
are similar system as Xe+Sn) . On the contrary, at higher 
energy (Fig. [TJc) and (d)), to reproduce experimental 
data, the weight of MDI should be increased. More in- 
formation could be extracted once more experimental 
data are provided to low down the statistical fluctua- 
tion. In addition, at energy near Fermi energy as shown 
in Fig. [TJa), it does deserve performing experiment to 
obtain Z-dependent stopping, which is important to un- 
derstand the minimal stopping at event level at Fermi 
energy as well as constrain the EOS. 

The atomic number dependence of stopping observ- 
able suggests that the heavier fragments have less stop- 
ping. This indicates that the heavier fragments are more 
transparent and keep a stronger entrance channel mem- 
ory even in the central nuclear collision. This supports 
the argument that the fragments are not formed in a glob- 
ally equilibrated environment (2~l| , which contradicts the 
assumption of equilibrium, required by ideal fluid hydro- 
dynamics. A self-consistent explanation is provided in 
[19( : if a heavier fragment can survive during the process 
of HIC, its constituent nuclcons should have suffered a 
less violent collision history. Basing on this considera- 
tion, the stopping hierarchy at fragment level would also 
reflect the N-N cross section. This explanation are also 
quantitatively demonstrated in earlier version of QMD 
[33 ] . Another explanation is due to corona effect (geo- 
metrical) [HI, even in very central collsion, those nucle- 
ons near the surface of projectile or target would suffer 
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less collision history than those inside. Thus, they would 
keep their original direction and compose the final heavy 
residues. 

In very central HIC, the incident energy dependence of 
stopping can be understood in a unified structure based 
on Pauli-blocking, mean field and N-N collision. At in- 
termediate energy, two-body collision between nuclcons 
always increases the dissipation of the HIC system and 
tends to favor more isotropy. While, Pauli-blocking sup- 
presses N-N collision at low energy, and mean field forms 
collective motion. This picture shows incident energy 
dependence. Below Fermi energy, the HIC system has 
enough time to reach an equilibration state because of the 
slow process of the reaction. Near Fermi energy, N-N col- 
lision is suppressed by Pauli-blocking greatly and nucle- 
ons tend to keep their original moving state and have an 
entrance effect. Meanwhile, mean field dominates HIC, 
and clusters are then favored. Above Fermi energy, two- 
body collision between nucleons comes to play the dom- 
inant role and reach their maximal effect at 300 - 400A 
MeV. While, at higher energy, the mean free path of nu- 
cleon increases because of the smaller N-N cross section, 
thus most of nucleons will pass through with each other. 

V. SUMMARY AND CONCLUSIONS 

In summary, we successfully describe, for the first time, 
the wide-range excitation function of nuclear stopping 
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